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Abstract 

The steady state existence problem for Kraichnan advected passive vector models is considered 
for isotropic and anisotropic initial values in arbitrary dimension. The models include the mag- 
netohydrodynamic (MHD) equations, linear pressure model (LPM) and linearized Navier-Stokes 
(LNS) equations. In addition to reproducing the previously known results for the MHD model, we 
obtain the values of the Kraichnan model roughness parameter ^ for which the LNS steady state 
exists. 
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I. INTRODUCTION 

This is a companion paper to a previous work by the present author wherein the 
phenomenon of anisotropic anomalous scahng was studied in the context of passive vector 
fields. The work was in part incomplete, as the main assumption was the existence of a 
steady state solution for the pair correlation function. We aim here to find exactly the 
preconditions under which this assumption is valid. Much of the technical material is from 
the above paper, to which we often refer for details. We study the stability of an equal time 
pair correlation function of a field u(t, x) determined by the equation 

iii — uAui + V ■ V-Uj — an ■ Vfj + VjP = 0, (1) 



where the vector field v(t, x) is determined by the Kraichnan model {2] and all vector quan- 
tities are divergence free. 

The model was introduced in jsj] as the most general linear passive vector model re- 
specting galilean invariance. The parameter a = —1,0 or 1 corresponds respectively to the 
linearized Navier-Stokes equations (abbreviated henceforth as LNS), the so called linear 
pressure model 3-8] (LPM) and the magnetohydrodynamic (MHD) equations jsl. loijlTl. In 
the context of the magnetohydrodynamic case, the inexistence of the steady state is known 
as the " dynamo effect" , where the dynamo refers to exponential growth in time of the pair 



correlation function (see e.g. |9l-lll| and references therein). This problem is by far the 



easiest of the three due to the vanishing of the nonlocal pressure effects. In j4| it was shown 
that for the linearized pressure model (corresponding to a = 0) the steady state always 
exists by showing that the semi-group involved with the time evolution is always positive. 
The analysis for the linearized Navier-Stokes case with a = — 1 is considerably more difficult 
than the above two cases since unlike in the MHD case, the nonlocal effects are present and 
contribute strongly to the dynamics, and because unlike in the LPM case, the semigroup is 
not always positive. 

The present goal is therefore to find the values of ^ for which the LNS steady state exists, 
where ^ is the roughness exponent of the Kraichnan velocity field w ~ r^. The method by 
which this is accomplished involves applying a Mellin transform on the eigenvalue equation 
for the two point correlation function, and by solving a resulting recursion relation. 
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II. THE MODEL 



We sketch here the derivation of the equation in Melhn transformed form and refer to 
the previous paper [l| for further details. All vector quantities in eq. ([T]) being divergence 
free results in an expression for the pressure, 

P = (1 - a) i-A)-' diV.djUi. (2) 

One may then rewrite the equation compactly as 

iii - uAui + Vijk (ujVk) = 0, (3) 

with an integro-differential operator 

T^ijk = Sijdk - a5ikdj + (a - l)(9i9j9fcA"\ (4) 

where is the inverse laplacian. The equal time pair correlation is defined as 

Gij{t, r) = {ui{t, X + r)uj{t, x)), (5) 

where the angular brackets denote an ensemble average with respect to the velocity field, 
which in turn is defined by the Kraichnan model as 

(i;,(t,r)t;,(0,0)) = 5{t)D,,{r) 



^iqr 



where we have defined the incompressibility tensor Pij{q) = Sij — qiqj, df^q := -4:4^ and 



(27r)° 

^ G [0, 2] is a parameter describing the spatial roughness of the flow. 
We note a subtle difference from ^ in that we define the constant 

The reason for this is that the velocity correlation and structure functions would otherwise 
diverge at ^ = and ^ = 2 as the mass cutoff is removed. This aspect of the Kraichnan 



model has been clearly discussed in 25|. The equation for the pair correlation function is 
then 

dtG.j - 21. AG,, - V,^,V,p^ {D^aG^p) = 0. (8) 
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The Fourier transform of the correlation function will then be decomposed in terms of 



hyperspherical tensor basis according to the prescription in 



m as 



(9) 



a,/ 



where the tensor basis components are 



Bffip) = \p\^-'d,d,<l>\p) 



,3,1 



Btfip) =\p\-\prd,+p,d,)¥{p) 



(10) 



?4,i 



[ B^Hp) = \p\-''^PiPj¥{p) 



and where is defined as := \p\^Y\p), where is the hyperspherical harmonic 

function (with the multi-index m = 0). It satisfies the properties 



^<^\p) =0 
p.V$'(p) =l^\p). 



(11) 



Note that we are concerned only with even parity and axial anisotropy. We now introduce the 
Mellin transform which will be used to transform the equation into a recurrence/differential 
equation. The method was (probably) first used in in the context of the hypergeometric 

20l | also has a useful section on the Mellin transform applied 



function. The textbook by Hille 



to differential equations. Appendix Elof the present work also contains helpful material, and 
also offers some insight into some limitations of the method. We define the Mellin transform 
of a Fourier transform of G (the anisotropy index / will usually be omitted) as 



ga{t, z) 



dw 



W^+'Gait.w) 



w 



(12) 



and the inversion formula 



Ga{t,r) = / dz\r\^A,ga{t,z) 
Jo- 



(13) 



with the definition Az = ^^^1^-*'"/^^^/?'' , which originates from the inversion of the fourier 
integral (with volume of the unit sphere Qd), and the subscript 0— was used to denote a 
contour from —zoo to +ioo passing z = from the left (see appendix|A]for details). We also 
often denote f{z) = Azg{z), which is just the ordinary Mellin transform. 



Applying the Fourier transform, dividing by p^, applying the Mellin transform and finally 
setting the cutoff parameter rriy to zero in eq. (see |]| for details), we obtain the complex 
recurrence/differential equation 

dtUt, z-2) + 2iyUt, z) - XUt, z - = 0, (14) 

with the definitions 

A = (a - 1) (rf + 1 + a(l - 0) rf<r(rf/2)Q ^ 

^ ^ 16sin(7re/2)r(^^ + 2)r(^^ + l) ^ ^ 

and 

™. _ 4er(^) r(rf/2 + /-«)r(rf/2-/3)r(a + /3-d/2) ^,,,^, 

^2a,2/3 r(l-e/2) T{a)Tif3)T{d + l~a-/3) ^ ^' ^ ^ 

where the matrix coefficients t"'^{z) are listed in the appendix of jj. We have also effectively 
set Z^o = 1 by redefining time and viscosity. Requiring the correlation function ([9]) to be 
divergence free, i.e. zero when contracted with pi, results in only two of the four coefficients 
Qa being independent. The resulting equation may then be written in the following form, 

dth{t, z + ^-2) + 2z/h(t, z + - f Al + A + B ■ x) h(t, z) = 0. (17) 



Here we have performed a translation z ^ z + ^, defined the vector quantity h = {gi,g2 
and the matrices by 



- d+^,d+z 




(18) 



A. Isotropic sector 

We will be mostly concerned with the isotropic case since much of the actual computations 
can be neatly performed all the way. For I = 0, only the tensors and i?^ are nonzero, 
and correspondingly in the tensor decomposition we only have the coefficients gi{z) and 
(jiiz) = —gi{z) (due to the divergence free condition). The equation f|T71) then becomes 
a scalar equation for gi{z) alone, hence we only need the (1,1) component of the matrix 
^Al + A + B ■ X j , which reads explicitly 
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A1 + A + B-X 



_ 2(a - l)(a^ - 1 - a - d)T (1 + ^/2) T (1 + d/2) 
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r(^) 



r(-./2)r(^) 
■2r(^)r(^) 



-P^(^) .^..JJx^\4Vn = (19) 



where 



= -(a 1)^(2-0 
+ (2 + e - 2) {{d - + (^(^ _ 1) + 2aO ^+ 

C (-rf - 1 + 2a(d + 1) - a^(l + 2d - + ^ _ (i^))) (20) 
The equation in the isotropic sector is then 

dtUt, z + ^-2) + 2z/^i(t, z + A|(^)^i(t, z) = 0. (21) 

III. THE METHOD 

We now consider the eigenvalue problem with gi{t, z) oc e~^^g{z), resulting in the equa- 
tion 

E-g{z + e - 2) - 2ug{z + + A?(^)^(^) = 0. (22) 



This is analogous to the Schrodinger method in |10l4l2l|. The steady state exists if one can 
show that the spectrum is nonnegative. However, for example in the magnetohydrodynamic 
case as in the above mentioned papers and in [9,], it was shown that there exists a critical 
value of the parameter ^ above which one has negative energies resulting in an exponential 
growth in time. This phenomenon is interpreted as the dynamo effect of magnetic fields. 
Previous studies on the dynamo problem have resorted to some approximative or numerical 
schemes to find the growth rate \E\ as a function of C,- Here we will settle for simply finding 
the values of ^ for which the energies are nonnegative, thus implying a steady state. This 
is done by studying the zero energy equation, i.e. setting E = 0. This has the advantage of 
providing us with an exact solution, up to a numerical solution of a transcendental equation. 
The argument used in the present work is closely analogous to the classic "node theorem" 



[see e.g. 



2l| and also appendix [A]), which can be roughly stated as follows. Suppose that 
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for some large enough value of E the corresponding solution /£;(r) is oscillating between 
positive and negative values ar large r and satisfies the boundary condition fsiL) = 0, for 
some large L (tending to infinity). If the spectrum is bounded from below, we know that by 
decreasing E the zeros of /L(r) will move to the right and satisfy the boundary condition 
for a discrete set of E. The value of E for which the smallest zero reaches L is then the 
ground state energy. Since we are interested in whether or not the ground state is positive 
or negative, we can instead study the zero energy equation and ask for which values of C, the 
solution crosses over from nonoscillating to oscillating. We know that the ^ = equation 
corresponds to diffusion, i.e. a nonoscillating zero energy ground state. As we increase ^, 
we may discover that the solution becomes oscillating at large scales, which would imply 
a negative energy ground state and therefore instability. The large scale behavior of the 
solutions is determined by the negative poles of the Mellin transform, so the problem is then 
reduced to finding these poles and determining if and when they become complex valued. 



IV. ISOTROPIC SECTOR OF LNS 



The stability problem in the linearized Navier-Stokes case is closely related to the laminar 



flow stability problem as described in §26 of 2j]. The equation is derived from the Navier 



Stokes equation by decomposing the velocity field into v(r) + u(t, r), where v is a stationary 
solution and u is a small perturbation, resulting in the equation 

Ui - uAui + V ■ Vui + u-Vvi + ViP = 0. (23) 

The question is then whether or not the laminar flow is stable under such perturbations. 
Here instead the field v is supposed to model a statistical steady state solution of the full 
Navier-Stokes turbulence, as prescribed by the Kraichnan model. We are therefore studying 
whether or not the Kraichnan model is an adequate steady state description of turbulence 
in terms of stability. For a = — 1 we now have 



. 4(rf + Or(l + e/2)r(l + rf/2) T{-z/2)T{^) 
^^'^ r(^) P-'^'>2T{^)T{'-^) 

where 



-{z + 0' + d{2 + iz- 3)z - (3 - 00) + 4;^'- (25) 
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The problem is obviously a more difficult one than in the MHD case due to the transcendental 
nature of the function A^. We can however expand it as an infinite product of zeros and 



poles according to the Weierstrass factorization theorem (see e.g. (23|). The function 
may then be rewritten as 



oo 



tJi - K) [z-K) 

where the + and — signs refer to zeros or poles that have respectively positive (or zero) 
or negative real parts, see Fig. ([2]). We also have the poles 6^ = 2k, = —d — ^ — 2k 
and s{z) is some unknown entire function on the complex plane and 6 is a. z -independent 
Weierstrass factor that enforces convergence of the infinite product. It can be derived by 
showing that asymptotically as z — ?■ ±oo, the poles and zeros of eq. behave respectively 
as ±2k + const, where the constant term depends on ^ and d. It may certainly be possible 
to derive bounds for s{z) by asymptotic analysis of eq. (El]), but since it can not contribute 
to the pole or zero structure of the solution, we refrain from doing so. We can also neglect 
the explicit form of the constant 6 for the same reason. The zero energy equation from (!22|) . 
rewritten here as 

9i^) =]^9iz + 0, (27) 

can then be solved by the same methods as in appendix El with the strip of analyticity 
requirement — ^ < 7le{z) < 0, resulting in 



g{z) =a^{z)^{z)i2 



V) 



r (f^) r (^2fc±i_2^ 



X TT e-^/^^ [ ) ^ (28) 



k>0 



r 



where "ifiz) satisfies the equation '^{z) = e^*(^^\I'(2; + whose solution is again an expo- 
nential of an entire function. The following subtlety concerning the above formula should 
be observed: we deliberately chose to use the form ~ l/r(l — x) instead of r(x), where 
the two are related by the Euler reflection formula r(— 2;/2)r(l + z/2) = —ir/ sin(7r2;/2). 
The reason is that only in this form the strip of analyticity remains pole free, as per the 
consistency requirement. For example using F (^—^^ the above result would introduce 
poles at z = a'l — ^n for positive integers n, that would eventually permeate the strip of ana- 
lyticity. However, in some cases as we increase ^, the poles will enter the strip of analyticity 



and render the solution incompatible with the strip of analyticity requirement. We will also 
demand that the solution converges to zero at imaginary infinities in order to justify the shift 
of integration contour (see appendix [X]). It seems quite difficult to deduce the asymptotic 
behavior from the above formula, but we can study it by an asymptotic expansion of the 
exact form of eq. f l24l) . The function A^(2;) behaves asymptotically as ~ at imaginary 
infinities, so the asymptotic version of the difference equation (127|) reads 

-g{z) = z'^-g{z^i) (29) 

up to some irrelevant constant term. The asymptotic solution is then g(^z) = a^{z)r{z)^ . 
Multiplication by r{~z/2) /T {{z + d)/2) in defining / introduces a pole at z = 0, which takes 
care of the boundary condition. Then we have asymptotically f{z) ~ (T^{z)e~i^y^^^^^^^~^^'^\ 
where z = x+iy. Fourier expansion of would then contain terms such as sin (^^nz^ , which 
would spoil integrability for < < 2, unless n = 0. Therefore cr^lz) has to be a constant. 
We see now that the poles of f{z) occur for non- negative integers n at z = ^ — 2 + 2k + C,n, 
z = 2{k — 1) (with A; > 0) and at 2; = — ^n, where only the latter affects the large scale 
behavior. We draw the important conclusion that the poles have no effect on the steady 
state existence problem. By looking at Fig. ([1]) we can see how the first few large scale 
poles behave in various dimensions. In two dimensions the leading pole aj" enters the 
strip at around ^ ~ 0.28 and is complex for all ^, which implies that there is no steady state 
at all in the isotropic sector, at least for ^ < 0.28 26|]. In three dimensions the poles aj~ 



and become complex at around ~ 0.345 and enter the strip of analyticity at around 
^ ^ 1.08. Similar behavior occurs with different in higher dimensions, until at d = 9 the 
poles stay real for all < ^ < 2. We have plotted the value of C,c in Fig- (E]) in dimensions 
2 ... 9 together with the magnetohydrodynamic case. The fact that in some cases for large 
enough ^ the strip of analyticity condition is violated could possibly mean that the steady 
state exists also for some large values of ^. We will however be content with studying the 
cases for which such a violation does not take place. 

One important lesson of the present section is that the " complexification" hypothesis of 
l|, |9| is indeed an indication of instability of the fiow, in the sense that the imaginary parts of 
the scaling exponents correspond to oscillations of the correlation function and are therefore 
responsible for the instability. The second lesson is that one only needs to be concerned with 
the negative zeros of A^(z) when considering the stability problem, since they become 
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FIG. 1: A plot of the a = —1 poles of the solution g{z) in the isotropic sector versus in various 
dimensions. The dashed curves denote the real parts of complex valued poles. There is an infinity 
of poles, but all the others not displayed here are real. Also, all the poles in dimensions d > 9 are 
real in < ^ < 2. 

the poles in the solution f{z). The positive zeros appear only as zeros in f{z). 
V. ANISOTROPIC SECTORS OF LNS 



The anisotropic sectors can be studied with the same methods as above. We will however 
refrain from performing the actual computations and simply extend what we have learned 
from the isotropic sector to the anisotropic case, namely that one simply needs to study the 
complexification of the negative poles of the solution fa{z) (we now have a matrix equation). 
The role of will now be taken by the determinant of the matrix in eq. ( !T7|) . instead of 
just the (1, 1) component. Since we already know that the flow is stable in the anisotropic 
sectors for a = 1 and a = (see e.g. % llSl), we concentrate only on the a = —1 case 
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FIG. 2: A plot of the critical value ^c, above which the flow becomes unstable, versus the space 
dimension d. The upper curve refers to the magnetohydrodynamic case with a = 1 and the lower 
to the linearized Navier-Stokes case with a = —1. 



in various dimensions. The anisotropic linearized Navier-Stokes exponents differ from the 
magnetohydrodynamic ones in that even if the leading exponent is real, the next to leading 
exponent may be complex valued, resulting in oscillating behavior at intermediate scales. 
These exponents however have no effect on the existence problem. This is because the 
boundary condition at R tending to infinity can only be satisfied by the leading oscillating 
exponent. We also need to make sure that the periodic function cr^{z) is again required 
to be a constant due to integrability, so that it will not cancel with any of the poles. This 
results from the fact that all matrix coefficients beside the (1,1) coefficient in eq. (lT9l) behave 



asymptotically as ~ irrespectively of / 27|, and therefore so does the determinant. Hence 
the same conclusions will be drawn as in the isotropic case, i.e. that is indeed a constant. 

We have plotted some leading negative exponents in various anisotropic sectors in two 
and three dimensions in fig. ([3]). In two dimensions the leading anisotropic exponent is 
real, except (strangely enough) for < 0.15. The anisotropic sectors are therefore quite 
stable in comparison to the completely unstable isotropic sector. In three dimensions the 
I = 2 anisotropic leading exponent becomes complex at ~ 0.937 and all the higher 
sectors have purely real leading exponents. The anisotropic sectors are therefore much more 
stable in comparison to the isotropic sector critical value ^c°^ ~ 0.345, somewhat similarly 
to the magnetohydrodynamic case. We also note that none of the poles lie inside the strip 
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of analyticity, so the results should hold for all ^. In dimensions d > A the anisotropic 
exponents are always real. 



VI. CONCLUSION 



The stability analysis of the passive vector models previously considered in a companion 
paper [l| was successfully completed. The critical value below which the steady state 
exists was found in all dimensions, although the possibility of a steady state for an even 
larger region could not be excluded in the isotropic sector. The reason for flow instability 
was shown to be caused by the complexification of the largest negative pole of the solution, 
corresponding to large scale behavior of the correlation function. It was observed that in two 
dimensions the linearized Navier-Stokes problem is not stable for any ^ > in the isotropic 
sector, but relatively stable in the anisotropic sectors. In three dimensions the isotropic 
sector was observed to be stable for ^ < 0.345, the I = 2 anisotropic sector for ^ < 0.937 and 
higher anisotropic sectors for all ^. In dimensions from four to eight, the isotropic sector is 
stable below the critical values plotted in Fig. (|2]) and the anisotropic sectors are stable for 
all ^. In dimensions d> 9, all sectors are stable for all ^. 
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Appendix A: Mellin transform and examples 



The Mellin transform has been used to solve differential equations previously in e.g. 



19 



20| . The purpose of the present appendix is to clarify some aspects of the Mellin 



transform method and to point out some of its limitations. 



1. Solving differential equations by Mellin transform 

Define the Mellin transform of a function /(r) with r > as 



oo 

f{z) = [ —W^fiw). (Al) 
J w 



We will be concerned with finite diffusivity /viscosity v in our equations, which amounts to 
/(O) = 1 (neglecting normalization), and power law or faster decay at infinity with some 
so far unknown exponent — The complex parameter z in the above formula is therefore 
restricted to —TZe{Q < TZe{z) < [29]. The inverse transform is then 

C+lOO 

/(r)= Jdzfizy (A2) 

C — lOO 

where —7le{() < c < 0. Because of the constant boundary condition at zero, f{z) must 
have a pole at z = 0. We will therefore take c = such that the contour will pass z = 
from the left, and denote this by 0— under the integration sign. We can now use the Mellin 
transform to define a differential/integral operator of order a as 

(5./)(r) = / dzA^iz)fizy-\ (A3) 
Jo- 

For example the derivative dr would correspond to a = 1 and Ai{z) = z. Consider now the 
equation 

{DJ){r) + Xfir) = (A4) 
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with the above mentioned boundary conditions. Expressing this with the help of the Melhn 
transforms yields 



where in the second term on the second line we have simply changed the integration variable 
z ^ z — a and shifted the contour from a — )■ on the third line (assuming sufficiently fast 
decay at imaginary infinities), and also denoted in the last sum the possible contribution 
of poles inside the strip —a < z < {]. We can solve the problem by solving the recursion 
equation 



but only if we can find a solution for which the strip of analyticity is pole free, i.e. that 
the sum of residues above vanishes, which also implies TZe{() > a. It turns out that in 
the present context there are situations where such solutions cannot be found, at least not 
without improving the present procedure. We will however be able to find a sufficiently 
large class of solutions for which this problem does not appear. 

2. Isotropic a -model equation for ^ = 2 

For ^ = 2 the problem becomes simple enough to be solved exactly even for nonzero 
energies. From eq. f[T^ we have now 




dz {A.(^)/» + A/(z - a)} r^-'^ + ^ r^-7^(/(^,)) = 



(A5) 



A,(z)/(^) + \f{z - a) = 



(A6) 



Kl{z) = {d- l)z^ + (4a + d{d - I)) z + 2a [d + a{d^ - 2)) 



(A7) 



and the equation (!22|) can now be written as 



{z-z+){z-z_)g{z) = 2ug{z + 2), 



(A8) 



where we have defined the roots 




(A9) 



z± = 



2{d-l) 
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with the discriminant 



De = -4(rf - 1)E + + d{d - 2) (rf^ - 8(d + 1)0^) . (AlO) 

-z/2) 



Employing the definition f{z) = Azg{z) oc (y'(z) and a translation z — > 2; — 2, we 



obtain the equation 



"^"^"^ ^ f{z) + ^fiz-2) = 0, (All) 



{z - z+ -2){z - z_ -2)' ^ ' 2u 
which should be compared to eq. flA6l) . A general solution to eq. flAlip can be written as 



where a2{z + 2) = (y2{z) is a so far arbitrary periodic function (with the subscript denoting 
the period). The solution (modulo the a2 -term) has poles at z = 2?t, and z = z^ — 2n 
for nonnegative integer n. The width of the strip of analyticity therefore has to be —2 < 
TZe{z) < 0, from which we conclude that cr2{z) must be an entire, periodic function and that 
the poles TZe{z±) < —2. At imaginary infinities the above solution behaves asymptotically 
as 

f{z) ~ (2z/)-"/V2(z)e-il^l2'^-^ (A13) 
Since (T2{z) is entire, we can expand it in Fourier series as cr2{z) = ^ CLk sin {kirz). Anything 

fc>0 

else than k = would however spoil the above asymptotic behavior, so we conclude that 
c"2(2) = const. Inverting the Mellin transform then yields 

G(r) - C, I j / (^) r (^) (rV2uf' 

G(r) - C\ JJz ^ ^ ^^^^2^ , (A14) 

which is what one would obtain e.g. in the MHD problem with a = 1 by a direct coordinate 



space solution joj. The MHD problem for ^ ^ 2 can also be easily solved, although we will 
not reproduce that calculation here. The ground state energy is then the value of E for 
which the discriminant De vanishes (crossover between oscillating and power law decay), 
i.e. 

d(d^- lOc/2 + 9rf + 16) ^ ^ ^ 
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which is negative for 3 < d < 8 implying a dynamo effect for ^ = 2. It is now tempting 
to use the above solution also for other values of a. However, for example for a = —1 and 
d — 3 we have 

= ± ^V-2E-15 (A16) 

which means that for energies E > —15/2, the poles z± are always inside the strip of 
analyticity. We must therefore conclude that in cases such as this, the method is not 
sufficient to determine the existence of a steady state. 
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FIG. 3: A figure showing the leading large scale poles in various anisotropic sectors in two and 
three dimensions. The dashed lines denote the real parts of complex valued poles. The thick curves 
denote two poles very close to each other. All the poles beyond what are shown here are real for 
< e < 2. 20 



